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Direct Flux, Inter-reflections, and Final 
Luminances in Rectangular Rooms 


By R. O. PHILLIPS, B.Arch., A.R.A.L.A. (Member). 


Summary 


A method is developed from first principles for finding the average 
illumination on any surface of a rectangular room due to the luminances 
of the other surfaces. This method depends on the Hisano “ room index,” 
k, and an associated function Pr, derived from the work of Yamauti, which 
is tabulated for a range of values of k. Illustrative examples show how 
it may be used to determine the direct flux required to produce a pre- 
determined luminance pattern, or to reverse this procedure. Some com- 
parisons are made between the results obtained by this method and by 
others, including direct measurement, and the agreement is shown to be 
generally satisfactory. 


(1) Introduction 


The report on the “ Design of the Visual Field ’(') opened for investigation a whole 
new range of ideas in lighting design. The stimulating paper by Waldram(2) carried 
some of these ideas a stage further. In one section of this paper he pointed out how 
the light released within an interior causes direct illumination of the surfaces, but 
that in addition inter-reflections build up a gross luminance pattern, i.e., the mean 
luminance of fairly large parts of the field. He later went on to show how the 
luminance pattern required could be specified by the concept of “ apparent brightness.” 
With the luminance of each surface predetermined, it is then possible to state the 
required total illumination on it. This may be considered in two parts: (a) The direct 
illumination; (b) The inter-reflected illumination due to the luminances of all other 
surfaces visible from it. If (b) can be found, by deducting it from the required total 
| the necessary value of (a) is known, which can then be treated as straightforward 

illuminating engineering. Waldram suggested that (b) could be found “often 
sufficiently accurately by inspection.” 

The later paper by Bellchambers and Ackerman(?) described an investigation into 
the same problem of finding the direct illumination required to give a predetermined 
total luminance. The authors pointed out that it is now possible, with the aid of the 
inter-reflection tables developed by Moon and Spencer, to obtain the final luminance 
of the room surfaces when the direct illumination is known, but that the reverse 
procedure is so complex as to be impractical. An attempt was made “to overcome 
these difficulties from a consideration of fundamental theoretical principles, but a 
method simple enough for practical application could not be found.” The authors 
therefore turned to an experimental determination of the ratio between total and direct 
illumination for a range of conditions. 

It is evident that design based on the predetermination of the luminance pattern, 
from which the required direct illumination is derived, could be valuable. However, 

* it seems to require a new approach to the problem of inter-reflections; neither evaluation 
“by inspection,” nor by experiment, necessarily restricted to a limited number of con- 
ditions, appears to provide a satisfactory solution. A theoretical treatment, sufficiently 
simple to be of practical value, should prove preferable. The present paper attempts 
to provide such a treatment which can be used either to find the illumination required 
to give a specified luminance pattern, or to reverse this procedure. 








a The author is on the staff of the New South Wales University of Technology, Sydney, Australia. “The 
+ manuscript of this paper was first received on August 29, 1955, and in final form October 17, 1955. 
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(2) Methods of Calculating Illumination 


Classical methods of calculation were based on the determination of the direct 
illumination at a point from a source of known intensity. The Harrison and Anderson 
““ Lumen Method ” determined the average illumination of a surface (the working plane) 
from sources of known flux, as modified by inter-reflection from other surfaces. The 
Moon and Spencer tables extended the range of the latter to include the determination 
of the luminance of room surfaces from sources of known flux. , y 

The inter-reflected component of the gross luminance pattern is determined by 7 
the average illumination of one surface from others of known luminance. There are | 
well-known methods for finding the illumination at a point from surfaces of given | 
luminance, and the average illumination could thus be found by considering a number 
of points distributed over the surface, but the procedure would be very laborious. 

A more direct method is to find the total flux on one surface from another; il 
this is divided by the area of the surface, the average illumination is at once deter- 
mined. The theory necessary for this determination seems first to have been worked 
out by Lambert(*), and has been further developed by Yamauti(5) and Moon(®). How- 
ever, since this work appears to have received little attention it will be briefly developed | 
here in Section 3.1. It will then be shown how, with the aid of the Hisano “ equivalent | 
square room,” a new function, Py, may be derived to simplify the work. 


(3) Theoretical Development 

(3-1) Total flux from one rectangular surface to another : 
Consider a luminous strip of width a, and of uniform luminance B*, extending 
vertically to infinity, as shown in Fig. 1. It is required to find the total flux from it to 
the horizontal surface LMNO, of width a and depth b. Taking O as the origin of 
co-ordinates on the horizontal plane, let P be the point (x, y). Then E, the illumination 
at P, is received within the angle OPL, so that :— 
OPL OPQ + QPL 


E=B. - =B 
27> 2r 


B x a-x 
=-——( arc tan- + arc tan 
2r Vv y 


The fiux received on an elementary area at P of dimensions Ay by Ax is E. Ay Ax. 
Therefore the total flux, ¢, on LMNO is given by: — 


ewe A a~—x 
= — arc tan +. arc tan dy. dx 
aT Jodo V y 


It can be shown that the evaluation of this integral gives :— 


d = Ba?.F (°) 
a 
b 


where F (*) denotes a “ flux function” of 


This function is defined for any variable vy by the equation+ 


* In this paper all surfaces are assumed to be uniformly diffusing, and all luminances are expressed © 
in foot-lamberts, or equivalent units. The inclusion in equation (2) of the often-troublesome factor #, then 
eliminates it from all subsequent calculations. 

+ By putting w = arc cot v, the numerator of the above expression may be transformed into the bracketed 
portion of @(w) given by Yamauti (5) and quoted by Moon (®) and Hisano (7) as:— ; 


Gfw) = bf w.cotw fin. sin. © — § cot®w. In. cos w} 


76 Trans. Ilum, Eng. Soc. (London). © 
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Tt 
If the luminous strip only extends down to the line RS, at a height h above the 
horizontal surface, the flux on the latter is given by :— 


[ Veh 
d = Bar 4 a -F(*) \ Lies Dh shhcone<<taxouecsde cee (3) 
L a a , 


Denoting surface RSLO as 2, and LMNO as 3, then ¢,,., the flux on 3 from 2, is given 
by the difference between equations (1) and (3) so that :— 


b h /pr+ pf! 
$3» = Ba2 f r( )+F( y-C —) \ Sib 16-4 Wise a sage ete Was aca tale aa (4) 
( a a a + 


It is to be noted that this expression applies equally well to the flux from LMNO to 
RSLO. This “reciprocity” of the total flux from one surface to another, for the 


» same uniform luminance, is a perfectly general property, as shown by Moon(®). 


If the luminous surface is parallel to the receiving surface as in Fig. 2, ¢,,, the 
flux from surface 1 to surface 3 is given by :— 


Methods for dealing with more complex cases are given by Moon. With equations (4) and 
(5) and a table of values for the “ flux function,” the total flux received on any surface 


7 ofa rectangular room from any of the other surfaces can be determined. 


xpressed © 


w, then 


racketed | 


London). 


(3.2) Flux in the square room 

The use of the above equations wouid still be rather laborious in most cases; 
fortunately they can be considerably simplified for a square room. Consider a square 
toom of side a and height h and let k equal the ratio between wall and ceiling areas 
so that: — 


area of one wall ah h (6 
. GE NIE LRT ee EN ere He ) 
area of ceiling a® a 
Vol. 21, No. 4, 1956. n 


XUM 





R: O. PHILLIPS 


and in such a room let P; equal the proportion of the total flux emitted by the ceiling i 
which reaches the walls. Let subscripts ,, ,, , refer to the ceiling, walls and floor [ 
respectively; then the luminance of the ceiling is B,, etc. The total flux emitted by the © 
ceiling is B,. a, and the flux from ceiling to all four walls is given by putting b = a in 
equation (4) and multiplying by four. So that :— 


Teeetotne a lt C3 a at aD 


Ba 


} : 
and as =k, then: P, = 4[ F() + F(k) — F(V1 4. h) |" 
a 


Hisano(’) showed by “ flux functions” that, to a reasonable degree of accuracy, 
the flux distribution in any rectangular room could be treated as that in the equivalent 
unit square room of height k. For a rectangular room of horizontal dimensions a by 0, 
and height h, k is given by :— 
average wall area total wall area 

floor area ~ 4+ floor area 


h(a+b) 

~ 2ab 

Moon and Spencer(®) state that equation (8), the most general expression for k, appears 
to apply to all rooms with flat ceilings and no re-entrant angles. 


ork = 


(3.3) The unit square room 

We are now able to find expressions for the average illuminations on all surfaces | 
of a room in terms of k, P,, and of B,, B,, and B,, the luminances of ceiling, walls, | 
and floor. 

As the ceiling is of unit area, the total flux from it equals B,. Therefore the flux © 
on the walls is B,. P,, and since their area is 42, E,, the average illumination on the | 
walls due to the ceiling is :— 

Ps 


E B.. 


21 — 
By symmetry, the flux from the floor to the walls produces an average illumination © 
of : — 
Py 
E., =B 
23 3 4k 
By the reciprocity principle, the flux from walls to ceiling is equal to the flux from 7 
ceiling to walls, for equal luminances. Therefore the average illumination of the — 
ceiling due to walls of luminance B, is: — q 
E 
The average illumination on the working plane at a distance h’ below the ceiling, 7 
due to inter-reflections, may be treated by considering this plane to be the floor of a” 
unit square room of height k’ where : — 


-~ 


rn 


By similar methods the illuminations of the other surfaces may be determined: _ 
the results are summarised for easy reference in Table | As will be shown bi e 


~ * This equation may be expressed directly in terms of k, as:— 
P, ee [4 V1-k2. arc cot V1 ~ k2 — 4k are cotk + kf in + k2)— Ink®(2 4 K))] 
us j 


OQ. 3k 
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DIRECT FLUX, INTER-REFLBCTIONS AND FINAL LUMINANCES 


Table x 


Notation and basic relationships between room dimensions, primary flux, 
and average illuminations and luminances of the room surfaces 


Room dimensions axb k Area of walls h(a +- b) 
~ 4 < Area of floor — 2ab s 
Height, floor. to ceiling h | 
} lide | 
Height, working plane to ceiling h’ a= h -R 
p Flux from unit square ceiling to walls of height k. | 
ats Total flux from ceiling el 
1 )—— , | 
= ]— [ 4V 1+k? arc cot. V1 +k? 4k. arc cot. k {2 In. (1 +k?)—In. k® (2 +a }] 
- 
-— —-—---— | 
Ceiling Walls Floor 
Area (sq. ft.) | Sy=a x 6. Sg=2h (a + 6). Sg=a X b. 
Reflection factor Pp P2 P3 
Primary (direct) flux ; 

(Im.) $'; ¢' > 3 | 
Primary illumination 
(Im. /ft.?) EF’, =’, S; b'5=9'9/So EB’; =9'5/S3 | 
Primary luminance 
(ft.-1.) B =p, Ey B’g= pz F's B’3=pg E's 
Secondary illmination ; | 
from: | 
Ceiling of luminance B, — E9;=B,. Pp/Ak E3,=By,. (1—Pp) | 
Walls of luminance B,  Ey2=Bg. Pp Eog= Bg (1—2P,/4k) | Egg=By.Pp 
Floor of luminance B, E43= Bg (1—P) Egg= Bg. Pp/Ak. | ao 
at 3 . ace as sae —| 
Total secondary illumi- | 
nation EF” ,=E 49+ E43 E"g==E9,+Eo9+ Eq, E’3=E3,+ Ege | 
Secondary luminance | B”,;=p,k", B",=pfok"s B’3=p3E"s 
Illumination on work- | 
ing plane due to ceiling 
and wall luminances  E,,,;=B, (1—P,’). Eug=By. Px’ 


examples, they permit the calculation of either the direct illuminations to produce 
given luminances, or the luminances produced by a known flux. 

It is necessary to be able to find the total direct flux (or average illumination) on 
each major surface. There is still a need for easy methods of achieving these results 
but Einhorn’s “ Sector Flux (9) or Potter and Russell’s “ Zonal Flux ”(!°) methods 
appear to offer possibilities. 

Table 2 lists values of P, for k up to 2, which will rarely be exceeded. The figures 
have been derived from a table given by Yamauti for his ¢, (w), by interpolation with 
second differences; it is possible that the last figure may err slightly in some cases. 


(4) Practical Application 
The equations developed in Section 3.2 permit the solution of the inter-reflection 
problem as suggested by Waldram, but in a more positive manner than by inspection. 
They also make adjustment between flux, luminance and reflection factors to give the 
desiréd result quite straightforward, as shown in the example which follows. The 
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question of maintenance factors, which are simply corrections applied to the design 
values for loss in efficiency, is not considered.. 


Example 

Suppose that a doctor’s waiting room measuring 12.5 ft. by 12.5 ft. and 12.9 ft. 
high is to have totally indirect lighting, giving at least 10 Im/ft? on a plane 2.5 ft. above 
the floor. Taking a value of 0.65 for the average reflection factor-of printed paper, 
the luminance of a magazine page will be about 6.5 ft.-lamberts. To avoid a change of 
adaptation when looking up at the walls, these will be designed for 6 ft.-lamberts. The 
floor is to be not less than one-quarter of this, requiring 1.5 ft.-lamberts. For efficiency, 
the ceiling reflection factor is selected as 0.8, while a restful wall colour tentatively 
chosen has a reflection factor of 0.4. Since the floor receives only indirect illumination, 
its reflection factor will be selected after the illumination on it is found. 








"Riise. 














area of walls 645 
k —_ ove = 1.033 
4 x area of floor 4x 156 ‘ 
p 
h’ 12.9-—2.5 i 
k’ R : x 1.033 = 0.833 
h ie 
From Table 2, P;, = 0.809 P,’ = 0.749 
P, P, 
= 0.196 1—2 = 0.608 
i ak 
To provide a minimum of about 10 Im/ft?, design will be based on an average of — 
12 Im/ft?. The illumination on the working plane, due to walls of luminance 6 ft.- © 
lamberts, is é 
Ey, = B,. P,’ = 6 x 0.749 = 4.5 Im/ft? 
Illumination to be provided by ceiling = 12 — 4.5 = 7.5 Im/ft? y 
. Ey, = B,. (1 — P,’) = 7.5 | 
a = 29.8 ft.-lamb ; 
i, 0251 8 ft.-lamberts. é 
A value of 30 ft.-lamberts will be used. The calculations can now be conveniently set F 
out as below: — i 
| | Ceiling | Walls | Floor | | 
oe - - | _ Y 
1 | Area 156 645 156 3 
2 Reflection factor 0.80 0.40 ? | 
3 | Required luminance | 30 6 1.5 Ff 
4 | Required final illumination, (3) + (2) | 378 | 15 
foal fa eae : | ” ie | 
5 | Illumination due to ceiling luminance 5.88 5.72 | 
6 | Illumination due to wall luminance | 4.85 3.64 4.85 | 
7 | Illumination due to floor luminance | 0.29 0.29 
= ianetiiainiingitnatoniaitetaiuens OP SEE NT ree | 
8 | Total indirect illumination, (5) + (6) + (7) 5.14 9.81 10.57 
9 | Required direct illumination, (4)—(8) | 32.36 5.19 | 
10 | Required direct flux (lumens), (1) « (9) | 5050 | 3280 | 











As the illumination on the floor is 10.57 Im/ft2, and its luminance is to be 
1.5 ft.-lamberts, the required reflection factor, p,, is 0.14. 
It is known that if four luminaires of a certain type are suspended 2 ft. from 


80 Trans. Hlun. Eng. Soc. (London), 
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design Table 2 
Values of the function P,. 
k. Py. diff. kh. Py. diff. ao cae .' oe 
0.01 0.0196 594 0.51 0.5914 g¢ 1.01 0.8028 9¢ 1.51 0.8905 |, 
2.9 ft. 02 0385 188 52 5979 6 1.02 8053 32 152 9916 1} 
above | «03S s«0569 188 53 6042 Bo 1.03 8078 32 1.53 3927 It 
vaper, #040748 178 54 © 6105 6 1.04 8103 3° 1.54 8938 1} 
al 05 0822 125 55 6166 8h 105 127 74 1.55 949 
The | 06 1091 15 56 6226 54 1.06 8151 94 1.56 8960 |, 
iency, | 07 1256 16° 57 6285-3 1.07 8175 38 1.57 971 Ii 
tively = & ww a 7 a ae ee 
, 0s = 59 6400 2 8221 33 g992 1) 
ation, 10 1730-134 60 6455 Be 1.10 8243. 32 1.60 9002 1° 
dt 1881 61-6509 1.11 8265 161 9012 
s oo 62 6563 34 1.12 8272? 162 9022 10 
= §«6&.- 63 6616. 3 1.13 8308 31 1.63 9032 |? 
14 2313 148 64 6667 2 1.14 8329 5} 1.64 9042 10 
15 2451 138 65 6718 34 115 8349 30 165 9052 '0 
16 2586 .~ 66 «6767 1.16 8369 1.66 9061 
v @S 67 «6815 48 1.17 8389 70 167 9070 
18 2847-15 68 6863 45 1.18 8408. 12 1.68 9079 3 
19 2974 127 69 ©6910 42, 1.19 8427 19 1.69 9088 2 
20 © 3098 524 70 © 6956 4° 1.20 8446 10 1.70 9097 8 
te 71 7001 44 121 8465 1g 71 9106, 
ge of § 0-22-3338 112 72 «7045 48 1.22 8483 j8 172 o114 4 
ft. § 23 3454 114 73 7088 43 8501 1 , : 
568 IIs 4 867131 43 124 6 asi9_ 18 1740 ogi 8 
b 25 © 3680 jhe 75 «(7173 4 1.25 8536 17 175 91398 
8 2 76 7214 4g 1.26 8553 7 1.76 9147 
P 27-3896 502 7 22 1.27 8570 {7 1791558 
; 28 4001 108 78 7293 30 1.28 8587 12 1.78 «91638 
29 4104 108 79 7332 32 1.29 8603. 16 1.79 «9171 8 
y 30 = 4205 10! 80 7370-38 1.30 861918 1.80 9179 § 
ae 31 4303 81 7407 . 1.31 8635 1.81 9186 
“5 2 wm = 827444 37 1.32 8650 {2 182 gia § 
33 4495 (O° 83 7480 3° 1.33 8665 {2 1.83 9201 7 
— 344588 - = 1.34 8680 12 1.84 9208 7 
ri 35 4679 85 7550 3° 1.35 8695 |> 185 9215 2 
r 36 4768 gn 86 7584 44 1.36 8710 1 1.86 9223 
; 374855 (8? 87 7617: -33 1.37 8725 |? 1.87 9230 7 
: a 38 494186 88 7650 33 1.38 8739 {4 1.88 9237} 
5 Fi 39 = $026 95 89 7682 35 1.39 8753 44 189 9244 5 
: 4 5108? 9 7714 3 1.40 8767 14 190 9251 2 
| 41-5188 91 7745 1.41 8780 1.91 9257 
2 | 42 5267 79 S. mo.w 2... a; 
SL 43 534578 93 7806 30 1.43 8806 13 193 9271 2 
b 44542178 94 = 7835 38 14419 13 194 9277 «§ 
| 45 5496 78 95 7864 30 145883213 195 9284 2 
j | n 46 556925 9% 7893 9g 1.46 8844 4 196 9290 
a 47 564172 97 7921 38 1.47 8857 }> 1.97 9296 6° 
| § 48 5711 98 7948 55 1.48 8869 15 198 9302 ¢ 
Ig 49-5780 99 7975 3 1.49 8881 15 1.99 9308 © 
> be 50584868 1.00 8002 3? 1.50 8893 1? 2.00 9314 
‘rom 
wdeuilt Vol. 21, No. 4, 1956 $1 
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the ceiling, 0.746 of their total effective output will reach the ceiling. Therefore :— 
: ; 5,050 
R i c= i ee 
equired outpu 0.746 
Of this total, the remaining 0.254 will reach the walls, so that :— 


Direct flux to walls = 6,770 x 0.254 = 1,720 lumens. 

This is roughly only one-half the required value. However, the required luminance 
of 6 ft.-lamberts on the walls can be maintained by adjusting their reflection factor. 
If this is done; we will still have the same indirect illumination of 9.81 Im/ft?. The 
average direct illumination. will now be:— 

— = 2.66 lm/ ft? 
giving a total of 12.47 Im/ft?. To obtain a luminance of 6, we require a reflection 
factor for p, of 0.48. 
A value of 0.5 will be adopted, which should result in a luminance of 6.24 ft.-lamberts. 
The final average illumination on the working plane will then be: — 
From the ceiling, E,,, = 30 x 0.251. = 7.53 
From the walls, E,,. = 6.24 x 0.749 = 4.67 
Total, E,, = 12.20 Im/ft2. 

This completes the design procedure. It can be seen that it is quite flexible, and 

that adjustments may be made without re-calculating from the beginning. 


6,770 lumens. 


(5) Inter-reflections from Primary Flux 


The equations of Table 1 can also aid in finding the illumination and luminances 
which result from a known primary flux. However, several methods are possible and 
some additional factors must be considered, as discussed below. 


(5.1) Method based on Table 1 


The data of Table 1 would permit us to trace reflected flux through stage after 
stage of the inter-reflection process. For example, with a direct illumination E’,, 
on the ceiling of reflection factor p,, its luminance B’, = p,E’.. 

Then the luminance B,,, of the floor, due to the ceiling of luminance B’,, is: — 

Bs, = pps E’, (1—Py) 
This value of B,, will in turn produce a luminance on the walls of : — 
’ _, 
Bog = pyCo03 EF’ (1—P) dpe (13) 

It is evident that this process could be continued indefinitely for each surface in 

turn, but that it would rapidly become quite unmanageable. 


(5.2) Summation by integrating sphere theory 

Alternatively we can turn to the theory of the integrating sphere, as has been 
done by Dresler(!!), to determine the reflected component in daylight. Suppose that 
flux ¢’ is released in an enclosure of surface area S, and uniform reflection factor p. 
After one reflection the added flux will be ’p, after two reflections 4’p2, and so 


on. Therefore the average illumination, E, after an infinite number of reflections will | 


be found by adding these terms and dividing throughout by S, giving: — 
ce? ,%p, Hp , 
ij Ss Ss Ss 

This is a geometric series, and since p is by nature positive and less than unity, the 
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sum of ali terms after and ‘including any given ‘term is’ found by dividing that term 
by 1—p, so that: — 


¢’ 
E = a EE i ade adhstceldS QRMGIS AS RRs Janu (5a) 
S(1—p) 
wr me Be Hit ct seal acti waste dee alae (15b) 
S S(1—p) 
-. Fy. oe, FP" 15. 
or. &.. = st S 4 MC iit £23 sek das oon tee ae eae (15c) 


To apply these equations to a room where the various surfaces have different reflection 
factors, it is necessary to substitute for p a value p, the mean reflection factor. 

Equation (13) shows that primary flux reflected in turn from the ceiling and the 
floor produces on the walls a luminance given by an expression containing the 
product p,p.9,. We can substitute for this (p)> where p is the mean reflection 
factor at each reflection. It is evident that in this case p is the cube root of p;peps, 
ie. it is the geometric mean of these quantities. 

It seems reasonable to assume that, after many reflections, the number of times 
any small element of flux has been reflected from each of the surfaces will be 
proportional to the relative areas of ceiling, walls, and floor, i.e., to 1, 4k, and | 
respectively. Then the weighted geometric mean is given by:— 


log p, + 4k. log p, + log ps 
Lie Ce su (16) 


Substituting the value of p so found for pin equation (15a) would therefore give 
at least a first approximation to the average illumination of all surfaces in the room. 
However, where the reflection factors of the various surfaces differ by ordinary 
amounts this method is too inexact. 

Improved results can be obtained by applying the summation process to the 
second term, as in equation (15b). It is convenient to adopt the convention of adding 
() to all photometric quantities due to primary flux, and (”) to those due to secondary 
(i... once reflected) flux. The term ¢’/S, which represents the primary illumination, 
is then treated separately for each surface, giving: — 


Sar Gis calng, BOS) hlvcoes iscawigs cosines fwd. 203, ccuneas (17a) 
foe the walle, Be PIG: ites etic ie cas asiicend dpnde Gil wake (17b) 
for Bie heor, = Bi, IB. as Hae HIN s OG TEE Sees (17c) 


The numerator of the second term of equation (15b),¢’p, gives the flux which has 
been reflected once in a room of uniform reflection factor. In the non-uniform room 
this is obtained by adding the first-reflected fluxes from each surface. Hence the 
average added illumination due to flux which has been reflected once or more is 
obtained by adding the separate reflected fluxes and dividing as in equation (15b), by 
(1 —p) and by the total surface-area. 

The final average illumination on any surface is then found by adding this value 
to the primary illumination on it. 

If this method is applied to the example of Section 4, adopting the values of 
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primary flux and reflection factors determined in the design procedure, a comparison 
of the results with the designed values is shown below :— 


Design value __ Diff. 


Final luminance of ceiling (ft.-L.) .......... 33.2 30 + 10.7% 
= " » Walls ae Fh aee hewitt 5.9 6.24 — 5.6% 
o * . floor ee ae ee 1.27 1.5 —15.3% 

Illumination on working plane (Im/ft?) ...... 12.45 12.20 « +2.0% 


The agreement is reasonably good, but in a case such as this, where most of the 
reflected flux comes from one surface (here the ceiling), and hence would mainly 
illuminate the other surfaces, the summation process spreads it evenly over all surfaces. 
This results in values too high for the brightest surface, and too low for the others. 
In many cases errors of this order might be unimportant, but accuracy can be 
improved by delaying the summation process, as shown below. 


(5-3) Combined use of Table 1 and summation 
In this case we consider : 
(a) The direct illumination on each surface. 
(b) The indirect illumination on each surface due to first reflection from each 
other surface, by Table 1. 
(c) The later inter-reflected illumination. 


These are equivalent to the three terms of equation (1 5c). 

The average primary illumination on each surface is first found, and when 
multiplied by its reflection factor, gives its primary luminance. 

The secondary illumination of each surface due to the primary luminance of 
each other surface can then be found from Table 1. These secondary illuminations 
are added for each surface, and the result multiplied by the reflection factor, giving 
the secondary luminances B”,, B”, and B”, of ceiling, walls and floor. If these values 
are multiplied by their areas and added, we obtain the total flux after two reflections. 
Then ;:", the average inter-reflected illumination due to flux which has been reflected 
two or more times, is found by dividing the result by the total surface area and by 
(1 — p), found from equation (16), giving : — 


BY, + 4k.B”, + B” 
(2 + 4k) (1 — p) 


For each surface, the value of -” is muJtiplied by the appropriate reflection 
factor, giving the added luminance due to all later inter-reflections. The three values 
of luminance for each surface (primary, secondary, and inter-reflected) are then added 
to give the final luminance. The average indirect illumination on the working plane is 
found from the final luminances of the ceiling and the walls. any direct illumination 
being simply added to the indirect values. 


E” = a ee ere ee ree te eet a (18) 


Example 


As an example of this method, the values of primary flux and reflection factors 
selected in Section 4 will be used, as a check, in the reverse procedure. 


Basic Data 


k = 1.033 k’ = 0.833 
P, = 0.809 P,/ = 0.749 
(1-P,) = 0.191 (1—P,’) = 0.251 
P,/4k = 0.196 p = 0.438 (from eqn. 16) 
(1—2P,/4k) — 0.608 (1-p) = 0.562 
84 Trans. Hlum. Eng. Soc. (London), 
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CALCULATIONS 


' ' 


Ceiling Walls Floor 

1 | Primary flux , 5050 1720 — 
2 . Area 156 645 156 
3 | Reflection factor 0.80 0.50 0.14 
4 | Primary illumination, (1) + (2) 32.4 2.67 —~ 
5 | Primary luminance, (3) » (4) 25.92 1.34 — 

| 6 | Secondary illumination due to ceiling — ° 5.08 4.95 

| 7 | Secondary illumination due to walls 1.08 0.81 1.08 
8 | Total secondary illumination, (6) + (7) 1.08 5.89 6.03 
9 | Secondary luminance, (3) x (8) = 0.86 2.95 0.85 
10 | Luminance due to F” (below), (3) « F” 3.22 2.02 0.57 
11 Final luminance, (5) + (9) + (10) 30.00 6.31 1.42 

| Designed Value 30 6.24 * ae 
Difference from designed value 0 + 11% 5.3 % 


B", + 4k.B", +- B"s 0.86 +. (4.132 x 2.95) + 0.85 


ee SS ee 6.132 x 0.562 
= 4.04 Ilm/ ft? 
Illumination on working plane due to ceijing = B,(1—P,’)=300.251= _ 7.53 
Illumination on working plane due to walls = B,P,' = 6.31 x 0.749 = 4.72 
Total = 12.25 1m/ft? 
Design Value = 12.201m/ft? 


The agreement between the pre-determined design values and those calculated by 
this method is about as close as could reasonably be expected. Therefore the method 
gives consistent results whether proceeding from luminance pattern to required primary 
flux distribution, or vice versa. However, this does not demonstrate whether the results 
achieved would agree with those obtained by measurement or by other methods of 
calculation. This question will be examined in Section 6. 


(5-4) Variation to allow for grading of wall luminance 


These methods all give average values for the luminances of the various surfaces; 
in most cases this will possibly be adequate. However, with lighting systems which are 
mainly indirect, there will be a definite grading of wall luminance downwards from top 
to bottom. This can quite easily be found by dividing the height of the walls into n 
horizontal strips. If the distance from the ceiling to the bottom of the first strip is k1, to 
the second k2, etc., then the flux from the ceiling to the first strip is B, P,, and to 
the first plus the second is B, P,,, etc. Therefore the flux from the ceiling to the second 
strip alone is B, (Py2—P,,) and so on. The illumination on each is found by dividing 
the above values by 4k/n, the area of each strip. The initial procedure is as before, 
until the average final luminance of all surfaces is found. The illumination on one 
wall is then approximately the average value due to the other walls, plus values due to 
the ceiling found as above, plus illumination due to the floor. The latter may be 
treated as an overal! average, or in strips as for the ceiling; unless the floor luminance 
is high it will usually make little difference which method is used. 
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Example ; 

To illustrate this method, a somewhat improbable arrangement of a high room, 
not square, with a luminous ceiling has been chosen. This provides a severe test of the 
method and can also be compared with measurement, as shown in Section 6. 

A room 13 ft. by 20 ft. on plan, and 20 ft. high, is lighted by a translucent glass 
ceiling, illuminated from above so as to give a primary luminance, B’,, of 100 ft.) 
lamberts. The reflection factors of the glass ceiling, the walls, and the floor are 0.70,} 
0.58, and 0.20, and k = 1.27. : 

The average final luminances of the surfaces are found, exactly as in the example” 
of Section 5.3, to be: — 


Ceiling 110.4  ft.-lamberts. 
5.84 ,, %» 


The graded illumination due to the ceiling is then found as described above; 
Consider the walls as divided into five horizontal strips, each 4 ft. deep. In terms of} 
k, each strip is k/5 deep, i.e., 0.254, and the area of each strip is 4k/5, ie., 1.016) 
The work can be set out as below. 


Depth Equiv. Vales of ici Illumination Final 
below k Pr 

ceiling (ft.) 
(1) (2) (3) (4) 


due to ceiling illumination 
4 of strip 
(*) * 7016 (Im /ft2) 


(Px for strip) 


0 00 0.0 
0.3723 53.7 
0.254 0.3723 
0.2178 
0.508 0.5901 
0.1321 
().762 0.7222 
0.0821 
1.016 0.8043 
0.0526 
20 1.27 0.8569 


Column (4) represents the proportion of the total flux from the ceiling reaching ” 
the strip in question. The resultant illumination is found by multiplying by the final f 
ceiling luminance (110.4) and dividing by the area of the strip. The values in the 7 
last column are found by adding the illumination due to the final luminances of the 7 
other walls and the floor, 13.2 Im/ft?. The final luminances would, of course, be” 
obtained by multiplying by the reflection factor; the results have been left in the 7 
form shown for later comparison with measured values. 4 


(6) Accuracy of the Results 


It seems probable that the large number of variables involved in the general 
solution of the inter-reflection problem would make any completely rigorous approach © 
either impossible or so complex as to be quite impractical. Certain compromises § 
and approximations therefore appear inevitable; the best that can be hoped for is” 
that they will not lead to errors-large enough to be important in practice. ? 

The major approximation adopted in the method here described is that of 
assuming that all illuminations may be averaged over the surfaces on which they are 
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incident (except in the method of 5.4). An attempt will be made to assess the errors 
involved by comparison with other data. 

The original Harrison and Anderson data were frankly empirical. Recently 
Potter and Russell(!°) have re-measured “ Utilisation Factors”* using more refined 
techniques, comparing the measured results with those obtained by the Moon and 
Spencer “ Interflectance ’* tables. They demonstrated that, in general, the measured 
values are higher than those calculated. In the discussion following this paper, Dr. 
Spencer pointed out that the integral equations on which the tables are based required, 
for their solution, the use of “kernels” in which exponential approximations were 
substituted for the true values. 

For example in a room with perfectly diffusing luminous ceiling and black walls 
and floor, the “ interflectance,” f, (the ratio between flux to working plane, and flux 
from ceiling), is taken as:— 

f, = e-?*’ (with the notation used herein) 
Dr. Spencer goes on to point out that exact equations can be derived for specific 
room shapes, and gives for a square room an expression which is equivalent to (1 — P;’) 
of this paper, as of course it should be. Another compromise adopted in the Moon 
and Spencer tables is to standardise on a working plane at one-quarter of the room 
height. so that k’ is always 0.75k. The degree of error involved in this approximation 
can be seen below :— 


k 0.0 0.2 0.4 0.6 0.8 1.0 1.2 

k’ 0.0 0.15 0.30 0.45 0.60 0.75 0.90 
(1—P,’) 1.0000 .7550 5795 4504 3545 .2827 .2286 
e~ 2k 1.0000 .7408 5488 4066 3012 .2231 1653 
error 0 A9pec. -53 pc. -9.7 pc. -I15 pc. -21 pc. -27.5 p.c. 


The example of Section 5.1 was deliberately chosen so as to fit exactly one set 
of the Potter and Russell measurements (namely that of the second column of figures in 
Table V). The measurements extended only to the determination of the “ Utilance.”’ 
In the example of Section 5.3, this is: — 


12.25 x 156 
Total Flux o....0 
This differs from the measured value of 0.300 by --6 per cent. The value 
calculated from the Moon and Spencer equations for the exact room dimensions and 
the true height of working plane is 0.259, a difference of — 13.6 per cent. Comparison 
with others of the Potter and Russell data gives similar results as shown below: — 


E,, x Floor Area 


Potter and Russell Table and Column No. 


V,1 V4 V6 1V,7 

Type of luminaire Indirect Indirect Indirect Direct 
| 1.20 0.825 0.70 0.53 

Utilance, ‘‘interflectance’’ method 0.210 0.339 0.411 0.849 
Utilance, method of Section 5.3 ().242 0.349 0.423 0.850 
Utilance, measured 0.251 0.378 0.438 0.885 


Unfortunately no measurements of the surface luminances of the experimental 
room are available. A comparison has therefore been made with the normal Moon 
and Spencer tables, using the tabulated data nearest to those of the actual room of 
Section 5.3. The results are as shown at the top of the next page. Agreement is 
surprisingly close. 

A word should be said about the determination of p as the weighted geometric 
mean, in place of the weighted arithmetic mean, as used by other workers. The 


* Both these terms are equivalent to what the C.I.E. has now recommended should be called “ Utilance.” 
the ratio between flux to working plane and flux from luminaire (as distinct from flux from lamps). 
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Average Luminances of | § 

Illumination, Ceiling Walls Floor 

Im /ft? ft.-1. ft.-l. ft.-l. 
Due to 5050 Im. on ceiling 10.3 30.0 4.25 1.01 
Due to 1720 Im. on walls 2.0 1.5 2.15 | 0.28 
Totals 2s |}. | 6.40 Ls A 

| | 

pe SOE ea ee eee ee Tee i 

Values from Section 5.3 12.25 30.0 | 6.31 1.42 


theoretical justification for this has been given. Some comparisons have also been 
made between published tables and calculations of the type shown above, and better 
agreement has been obtained with the geometric mean. 

As one part of his extensive experimental investigation of daylight, Pleijel(!?) 
measured the illumination on the walls of courtyards. A model court was used in 
which the reflection factors of the surfaces could be changed, the “sky” consisting | 
of a sheet of flashed opal glass, illuminated from above and resting on the tops of the 
walls. This is precisely equivalent to a room lighted by a luminous ceiling, and one 
set of measurements has been chosen as the example of Section 5.4. 

A series of measurements was taken of the illumination at points on the vertical 
centre line of one long wall, first with all surfaces matt black, and then with the reflec. © 
tion factors of the ground 0.2, and of the walls 0.2, 0.35, 0.58, and 0.84. Calculations ' 
have been shown for the value of 0.58 for the walls; illuminations have also been 
calculated for values of 0.35 and 0.84. The results are shown in Fig. 3a. Since 
Pleijel expressed his measurements as daylight factors, i.e., as percentages of the 
illumination from a complete hemisphere of sky with the luminance of the actual sky, 
the calculated data have been plotted in the same form. The full-line curves are 
measured; the broken curves were formed by plotting at the centre of each strip of 
wall the calculated average illumination of the strip and joining the points so found. 
It can be seen that the agreement is reasonable, but that the shapes of the calculated 
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Fig. 3. Hluminations on the walls of a courtyard 20m x 13m x 20m high. The 
full lines are Plejel’s measured Vilnes, the broken lines were calculated. 
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curves are incorrect. In Fig. 3b, the results are replotted showing only the reflected 
component of the total illumination. This indicates the degree of error introduced by 


|} considering only the average illumination contributed by walls and floor. 


(7) Conclusion 


A method has been developed from first principles by which the direct flux, the 
inter-reflected flux, and the resultant illuminations and luminances on room surfaces 
can be related. This method pre-supposes that the direct flux to each surface can be 
determined; the data here required may not as yet always be available in readily 
usable form, but work on this aspect of design is proceeding in several countries. 
However, when this information is available, the remaining calculations require nothing 
more elaborate than the table of values of P, presented here and a slide-rule, as used 
for the illustrative examples. 

Several examples have shown that the method appears to give values which are 
self-consistent, whether working from final luminances to primary flux or vice versa. 
Comparison with values obtained by measurement and the Moon and Spencer 
“interflectance ” tables shows generally satisfactory agreement. Many more such 
comparisons would be required before it would be possible to say that one method 
of calculation was more precise than the other. When the designer has the information 
he requires in the form of published tab‘és there is thus no reason to go beyond them. 

However, the methods here presented do appear to have the following 
advantages : — 

(1) They are more suitable for use in design, that is in proceeding from the luminance 
pattern to the required flux distribution, as illustrated in Section 4. 

(2) For dealing with special cases, not covered by the usual tables, they are easier to 
apply. The only complex mathematical expression required is that for P,. This is 
related purely to the geometry of the room, and can be tabulated once and for all, 
as none of the other variables enter into it. The “ interflectance * method also contains 
complex equations into which, however, all three values of surface reflection factors 
enter, so that a change of one of them requires recalculation almost from the beginning. 

It is the author’s hope that other workers will test the methods of calculation 
described here against experimental data, to secure further evidence as to whethes 
the calculated and observed results agree as well as they appear to do in the admittedly 
limited number of cases here presented. 
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Annual General Meeting 


The annual general meeting of the 
Iiluminating Engineering Society was 
held at the Royal Hall, Harrogate, at 
2.30 p.m. on Friday, May 11, 1956. The 
president, Mr. A. G. Higgins, was in ihe 
chair. 

In presenting the report of the council 
and the accounts for the year ending 
December 31, 1955, the president drew 
attention to the increase in membership 
of 84 during the year. He emphasised 
the need for new members in ali classes 
and said that there were many organisa- 
tions in the industry who should be 
encouraged to become sustaining mem- 
bers. The president also referred to the 
activities of the centres during the year 
and to the new edition of the IES. 
Code. The hon. treasurer explained 
that the excess of income over expen- 
diture during the year was due to 
fortuitous circumstances and he drew 
attention to rising costs. 

The adoption of the report and 
accounts was proposed by Mr. J. S. Mc- 
Culloch, who also proposed a vote of 
thanks to the president, officers and 
members of council for their conduct of 
the society's affairs during the year. 
This was seconded by Mr. A. G. Smith. 
The report and accounts were declared 
approved. 

The president 


also said that the 


council had given considerable thought 
to a revised constitution for the society 
and it was hoped to bring this to the 
attention of members in 
future. 

Dealing with the election of officers 


the near 





ANNUAL GENERAL MEETING 


and members of council for the 1956-57 ; 


session the president said that no other 
nominations having been received the 
following were declared elected :— 
President: Dr. W. E. Harper. 
Vice-Presidents: Dr. R. G. Hopkinson, 
Mr. J. S. McCulloch, Mr. E. B. 


Sawyer, Mr. C. C. Smith, Mr. W. T. F. [ 


Souter. 
Hon. Treasurer: Mr. W. Robinson. 
Hon. Secretary: Mr. J. G. Holmes. 
Hon. Editor of 

Transactions: Mr. C. R. Bicknell. 

Members of Council: Mr. A. D. S. 

Atkinson, Dr. H. H. Ballin, Mr. D. E. 

Beard, Mr. H. G. Campbell, Mr. J. B. 

Harris, Mr. A. Pott, Mr. G. C. Small, 

Mr. J. F. Stanley. 

After the remuneration of the 
auditors for the current year had been 
agreed the meeting was declared closed. 


Annual Subscriptions 


At an extraordinary general meeting 
of the society, held at the Royal Hall, 
Harrogate, at 2.45 p.m. on Friday, 
May I1, 1956, resolutions amending the 
by-laws in respect of annual subscrip- 
tions and the class of overseas member 
were considered and approved. 

As from December 31, 1956, the class 
of overseas member will be discon- 
tinued and the following new rates of 


annual subscription will come _ into 
effect :— ; 
t S. 
i ais hin ete-de-k stale se Bd 4 10 
Corporate Members (except 
| re are oer P 3 10 
Students (ages 16-24) .......... 1 0% 


Sustaining Members (minimum) 7 7 


Trans. Illum. Ene. Soc. (London), 




















| 
é 















eeting 

Hall, 
riday, 
ng the 
scrip- 
ember 


> class 
iscon- 
tes of 

into 


cE «& 
4 10 
7 
10 | 


j 





0 
0 
- 


sm Ww 


London), 








Printed by 
Arcus Press, Lrp., 
Temple-avenue and 
Tudor - street, 
London, E.C.4, 


England. 














